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Abstract 

Traditional studies of combinatorial auctions often only consider linear constraints. The rise of smart 
grid presents a new class of auctions, characterized by quadratic constraints. This paper studies the 
complex-demand knapsack problem, in which the demands are complex valued and the capacity of sup¬ 
plies is described by the magnitude of total complex-valued demand. This naturally captures the power 
constraints in alternating current (AC) electric systems. In this paper, we provide a more complete study 
and generalize the problem to the multi-minded version, beyond the previously known ^-approximation 
algorithm for only a subclass of the problem. More precisely, we give a truthful PTAS for the case 
(j) (z [0, ^ — 5], and a truthful FPTAS, which/wWy optimizes the objective function but violates the ca¬ 
pacity constraint by at most (1 + e), for the case S (^, tt — (5], where (j) is the maximum argument of 
any complex-valued demand and e, <5 > 0 are arbitrarily small constants. We complement these results 
by showing that, unless P=NP, neither a PTAS for the case cj) G — 5] nor any bi-criteria approxima¬ 
tion algorithm with polynomial guarantees for the case when (jj is arbitrarily close to tt (that is, when S is 
arbitrarily close to 0) can exist. 


1 Introduction 


Traditionally, many practical auction problems are combinatorial in nature, requiring carefully designed time- 
efficient approximation algorithms. Although there have been decades of research in approximating combi¬ 
natorial auction problems, traditional studies of combinatorial auctions often only consider linear constraints. 
Namely, the demands for certain goods are limited by the respective supplies, described by certain linear 
constraints. 

Recently, the rise of smart grid presents a new class of auction problems. In alternating current (AC) elec¬ 
tric systems 1101, the power is determined by time-varying voltage and current, which gives rise to two types 
of power demands: (1) active power (that can be consumed by resistors at the loads) and, (2) reactive power 


‘This paper appears in ACM Transactions on Economics and Computation, Vol. 5, No. 1, Article 7, October, 2016. DOI: 
http://dx.doi.org/10.1145/2955089 Extended abstracts containing some of the results have been presented in the International Con¬ 
ference on Autonomous Agents and Multi-Agent Systems (AAMAS 2014) and a workshop at the International Conference on 
Computer Communication and Networks (ICCCN 2014) |16[. 
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(that continuously bounces back and forth between the power sources and loads). The combination of active 
power and reactive power is known as apparent power. The ratio between active power and apparent power 
is known as power factor. In practice, the physical capacity of power generation and transmission is often 
expressed by apparent power. Electric appliances and instruments with capacitive or inductive components 
have non-zero reactive power. However, most electric appliances and instruments are subject to regulations 
to limit their power factors Q. It is vital to ensure that the total power usage is within the apparent power 
constraint, given the maximum power factor of power demands. 


In the common literature of electric power systems IM, apparent power is represented by a complex 
number, wherein the real part represents the active power and the imaginary part represents the reactive 
power. Hence, it is often necessary to use a quadratic constraint, namely the magnitude of complex numbers, 
to describe the system capacity. The power factor is related to the phase angle between active power and 


reactive power. Yu and Chau introduced the complex-demand knapsack problem (CKP) to model a one- 
shot auction for combinatorial AC electric power allocation, which is a quadratic programming variant of the 
classical knapsack problem. 

Furthermore, future smart grids will be automated by agents representing individual users. Hence, one 
might expect these agents to be self-interested and may untruthfully report their valuations or demands. This 
motivates us to consider truthful (aka. incentive-compatible) approximation mechanisms, in which it is in the 


best interest of the agents to report their true parameters. In |25] a monotone ^-approximation algorithm that 


induces a deterministic truthful mechanism was devised for the complex-demand knapsack problem, which, 
however, assumes that all complex-valued demands lie in the positive quadrant. 

In this paper, we provide a complete study and generalize the complex-demand knapsack problem to the 
multi-minded version, beyond the previously known ^-approximation algorithm. More precisely, we consider 
the problem under the framework of (bi-criteria) (a, /3)-approximation algorithms, which compute a feasible 
solution with objective function within a factor of a of optimal, but may violate the capacity constraint by a 
factor of at most /3. We give a (deterministic) truthful (1 — e, 1)-approximation algorithm for the case </> G 
[0, (5], and a truthful (1, l-|-e)-approximation for the case </> G (|, 7r—()], where is the maximum argument 

of any complex-valued demand and e, <5 > 0 are arbitrarily small constants. Moreover, the running time in 
the latter case is polynomial in n and - (this may be thought of as an FPTAS with resource augmentation; 
see, e.g., 111 ^). We complement these results by showing that, unless P=NP, neither a PTAS can exist 
for the latter case nor any bi-criteria approximation algorithm with exponential guarantees for the case when 
(j) is arbitrarily close to vr. Note that the difficulty when f £ (|, tt] is mainly due to the fact that demands 
are allowed to have both positive and negative real parts, which can cancel each other; this allows an optimal 
solution to pack much larger set of demands, within the available capacity, than any polynomial time algorithm 


can detect. We remark also that |24 251 show no FPTAS exists for the case f £ [0, | — 5]. Therefore, our 


results completely settle the open questions in |251. 


1.1 Contribution 

In Table[T] we briefly list the inapproximability and the best known truthful mechanisms for the m-dimensional 
knapsack problem (mDKP), m > 2, along with our results for three classes of CKP (namely, demands with 
maximum argument f £ [0, |], (^ G [0, vr — <5], and <f> > tt — S' where <5 is polynomially small in n, while S' 
is exponentially small n). 

2 Related Work 


Finear combinatorial auctions can be formulated as variants of the classical knapsack problem |[^[^ 141. 
Notably, these include the one-dimensional knapsack problem (IDKP) where a single item has multiple 
copies, and its multi-dimensional generalization, the m-dimensional knapsack problem (mDKP). There is an 
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CKP[0, f] 

CKP[f -r5,7r-(5] 

CKP[7r-5',7r] 

mDKP 

Inapproximability 

No FPTAS |24l|25l 

No 

(a, l)-approx 
(Sec.gJ 

Bi-criteria In¬ 
approximable 
(Sec.|4]l 

No FPTAS 
(see, e.g., (T^) 

Truthful Mechanism 

PTAS (Sec.[531> 

(1,1-r e)-FPTAS 
(Sec.|^ 

(1,1 e)-FPTAS 
(Sec.|g 

None 

PTAS (8) 
Bi-criteria 

FPTAS 


Table 1: A summary of results 


FPTAS for IDKP (see, e.g., |[gi). 

In mechanism design setting, where each customer may untruthfully report her valuation and demand, it 
is desirable to design truthful or incentive-compatible approximation mechanisms, in which it is in the best 
interest of each customer to reveal her true valuation and demand Q. In the so-called single-minded case, a 


monotone procedure can guarantee incentive compatibility |211. While the straightforward FPTAS for IDKP 
is not monotone, since the scaling factor involves the maximum item value, 0 gave a monotone FPTAS, 
by performing the same procedure with a series of different scaling factors irrelevant to the item values and 
taking the best solution out of them. Hence, IDKP admits a truthful FPTAS. We remark that monotonicity 
may be not enough for the incentive compatibility in the general setting. More recently, a truthful PTAS, 
based on another approach using dynamic programming and the notion of the maximal-in-range mechanism, 
was given in Q for the multi-minded case. We will use the maximal-in-range approach in this paper. 

As for mDKP with m > 2, a PTAS is given in 0 based on the integer programming formulation, but 
it is not evident to see whether it is monotone. On the other hand, 2DKP is already inapproximable by 
an FPTAS unless P = NP, by a reduction from EQUIPARTITION |14|. Very recently, | fT8| gave a truthful 
FPTAS with (1 -|- e)-violation for multi-unit combinatorial auctions with a constant number of distinct goods 
(including mDKP), and its generalization to the multi-minded version, when m is fixed. Their technique is 
based on applying the VCG-mechanism to a rounded problem. Based on the PTAS for the m-minded multi¬ 
unit auctions developed in Q, they also obtained a truthful PTAS for m-minded multi-unit combinatorial 
auctions with a constant number of distinct goods. Intuitively, a valuation function is m-minded if it is 
completely determined by the values on m different choices; for simplicity we call this type of valuation 
multi-minded in the rest of the paper. 

In contrast, truthful non-linear combinatorial auctions were explored to a little extent. Yu and Chau 
p5) introduced the complex-demand knapsack problem, which models auctions with a convex quadratic 
constraint. An earlier paper p4| also introduced the same problem without considering truthfulness by a 
different name called 2-weighted knapsack problem. One can regard, the complex-demand knapsack problem 
with strategic considerations as an auction design problem, where users bid on complex-valued items, and 
a feasible solution allocates one item to each user such that the total magnitude of allocated items is below 
a certain threshold. Even though some of the existing techniques can deal with combinatorial auctions with 


convex non-linear relaxations (see, e.g., |19|), those techniques require bounded integrality gap and yield 
randomized truthful-in-expectations mechanisms. 


3 Problem Definitions and Notations 


In this section we formally define fhe complex-demand knapsack problem. We presenf firsf fhe non-sfrafegic 
version of fhe problem where we assume all paramefers are known beforehand. Then we describe fhe sfrafegic 
version where each user k declares his/her valuafion funcfion defined over a sef of declared demands. In fhe 
laffer selling, we consider fhe case where users could lie aboul Iheir valuation functions and demand sels in 


order lo opfimize fheir ufilily funclions (see Sec. 3.5 for a formal definilion of ulilify). Towards fhe end of 
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this section, we present an application of the complex-demand knapsck problem to power allocation in (AC) 
alternating current electric systems. 


3.1 Complex-demand Knapsack Problem (non-strategic version) 


We adopt the notations from p5| . Our study concerns power allocation under a capacity constraint on the 
magnitude of the total satisfiable demand (i.e., apparent power). Throughout this paper, we sometimes denote 
= Re(z^) as the real part and = Im(z^) as the imaginary part of a given complex number v. We also 
interchangeably denote a complex number by a 2D-vector as well as a point in the complex plane. \v\ denotes 
the magnitude of v. 

We define the non-strategic version of the complex-demand knapsack problem (CKP) with a set [n] = 
{1,..., n} of users as follows: 


(CKP) 


subject to 


max 

xe{o,i}" 


VkXk 

kG[n] 


^ ^ dkXk 

k£jV 


< c. 


( 1 ) 

( 2 ) 


where dk = + id\ G C is the complex-valued demand of power for the /c-th user, C G M+ is a real-valued 

capacity of total satisfiable demand in apparent power, and Vk G M+ is the valuation of user k if her demand dk 
is satisfied (i.e., Xk = 1). Evidently, CKP is also NP-complete, because the classical 1-dimensional knapsack 
problem (IDKP) is a special case. 

We define a class of sub-problems for CKP, by restricting the maximum phase angle (i.e., the argument) 
of any demand. In particular, we will write CKP[(^i, (/) 2 ] for the restriction of problem CKP subject to 

(pi < maxkeAf Si,rg{dk) < 4>2, where aig{dk) G [0, vr]. We remark that in the realistic settings of power 

R dP- 

systems, the active power demand is positive (i.e., df> 0), but the power factor (defined by |^) is bounded 
by a certain threshold, which is equivalent to restricting the argument of complex-valued demands. 

From the computational point of view, we will need to specify how the inputs are described. Through¬ 
out the paper we will assume that each of the demands is given by their real and imaginary components, 
represented as rational numbers. 


3.2 Non-single-minded Complex Knapsack Problem (strategic version) 


In this paper, we extend the single-minded CKP to general non-single-minded version, and then we apply 
the well-known VCG-mechanism, or equivalently the framework of maximal-in-range mechanisms |22|. The 
non-single-minded version is defined as follows. By a slight abuse of notation, we denote Vk{-) as a valuation 
function for the non-single-minded setting. As above we assume a set M of n users: user k has a valuation 
function Vk{-) ■ T) —]R_|_ over a (possibly infinite) set of demands V C C. We assume that 0 G V, Ufc(O) = 0 
for all k G J\f, and w.l.o.g., \d\ < C for all d gT>. We further assume that each Vk{-)is monotone with respect 
to a partial order “A” defined on the elements of C as follows: for d, / G C, d A / if and only if 

\d^\ > \n, M'l > l/'l,sgn(d^) = sgn(/i^),sgn(di) = sgn{f). 


(See Fig. [T]for pictorial illustration.) We assume 0 A d for all d G V. Then for all k G J\f, the monotonicity 
of nfc(-) means that Vk{d) > Vk{f) whenever d'^ f. 

The non-single-minded problem can be described by the following program (in the variables d^): 


(NsmCKP) max Vk{dk) 
keJV 


s.t. 


Y^dk <C 

keJV 

dk GV for all k G M, 


(3) 

(4) 

(5) 


4 







Figure 1: A pictorial illustration for the partial order “A”: di A d 2 and d^ < d^. 


where | J2kej\f^k\ = y (X^feeA/"+ (Ylk&M particular interest is the multi-minded version of 

the problem (MULTiCKP), defined as follows. Each user A: G AA is interested only in a polynomial-size subset 
of demands Dk C V and declares her valuation only over this set (that is, is a set of 2-dimensional vectors 
of size |i2fc| = poly(n)). Note that the multi-minded problem can be modeled in the form (NsmCKP) by 
assuming w.l.o.g. that 0 G Dk, for each user k G M, and defining the valuation function Vk{-) ■ D —M+ as 
follows: 

Vkid) = max {ufc(dfc) : dk <d}. (6) 


We shall assume that the demand set of each user lies completely in one of the quadrants, namely, either 

> 0 for all d G Dk, or < 0 for all d G Dk- This assumption is needed for the results in Sec.[^ as we 
will see, the problem is split into two independent 2-DKP problems based on the demands’ quadrants, then 
each user k is allocated by one demand in the quadrant where all her demands Dk lie. (As we will see in 
Sec. 3.6 d^ > 0 corresponds to an inductive load, while d^ corresponds to a capacitive load.) Note that the 
single-minded version (which is CKP) is special case, where \Dk\ = 1 for all k. When we consider strategic 
users, we will assume that the users can lie about their demand sets and/or valuation functions (as long as the 
demand set of each user lies completely in one of the quadrants). 

We will write MULTlCKP[i;()i, cj) 2 ] for the restriction of the problem subject to cj)i < cj) < cj )2 for all d G P 
where cj) — maxdgx> arg(d) (and as before we assume arg(d) > 0). 


3.3 Non-single-minded Multidimensional Knapsack Problem 

To design truthful mechanisms for NsmCKP, it will be useful to consider the non-single-minded multidimen¬ 
sional knapsack problen:|^(NSM-mDKP) defined as follows, where we assume more generally fhat V C 
and a capacity vector c G is given. As before, a valuafion funcfion for each user k is given by Q. An 
allocation is given by an assignment of a demand dk = (d]., d™) G 22 for each user k, so as to satisfy the m- 

dimensional capacity constraint J2k&Af < c. The objective is to find an allocation d = (di,..., d^) G 
so as to maximize the sum of the valuations Ylk&jy '>^k{dk)- The problem can be described by the following 
program: 


(NSM-mDKP) max Vk{dk) 

(V) 

keJV 

s.t. dk < c 

(8) 

keJV 

dk €T> for all k ^ M. 

(9) 


'Sometimes, this is also called the multiple-choice knapsack problem. 
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Similarly, we consider the multi-minded version of the problem (MULTl-mDKP): each user fc G AA is 
interested only in a polynomial-size subset of demands ^ V and declares her valuation only over this set. 
The multi-minded problem can be modeled in the form NSM-mDKP by assuming w.l.o.g. that 0 G D^, for 
each user k G J\f, and defining the valuation function Vk{-) : T) —)• M+ as Vk{d) = max.d^^Dt,{vk{dk) ■ dk ^ 
d}, where is a component-wise partial order. 

It is worth noting that NSM-mDKP is similar to multi-unit combinatorial auctions (CA) with m distinct 
goods; the difference is that in the latter problem the set V is restricted to be integral, whereas we do not 
assume this restriction in NSM-mDKP. 

3.4 Approximation Algorithms 

We present an explicit definition of approximation algorithms for our problem. A feasible allocation satisfying 
Q is represented by a vector d = (di,..., dn) G P”. When no demand is allocated to user k, we assume 
dk = 0. We write r;(d) = optimal allocation of NsmCKP (or MULTiCKP) 

and Opt = r;(d*) be the corresponding total valuation. We are interested in polynomial time algorithms 
that output an allocation that is within a factor a of the optimum total valuation, but may violate the capacity 
constraint by at most a factor of /?: 

Definition 3.1. For a G (0, 1] and (3 > 1, a bi-criteria [a, (3)-approximation to NsmCKP is an allocation 
{dk)k & T3>'^ satisfying 


such that 




</3-C 


vkidk) > a • Opt. 

k£jV 


( 10 ) 

( 11 ) 


Similarly we define an (a, /3)-approximation to MULTiCKP. 

In particular, a polynomial-time approximation scheme (PTAS) is a (1 — e, 1)-approximation algorithm for 
any e > 0. The running time of a PTAS is polynomial in the input size for every fixed e, buf fhe exponenf of 
fhe polynomial may depend on 1/e. An even sfronger nofion is Si fully polynomial-time approximation scheme 
(FPTAS), which requires fhe running lime fo be polynomial in bofh inpul size and 1/e. In Ihis paper, we are 
inferesled in an FPTAS in fhe resource augmentation model, which is a (1,1 -|- e)-approximalion algorilhm 
for any e > 0, wilh fhe running lime being polynomial in fhe inpuf size and 1 /e. We will refer fo Ihis as a 
(1,1 -h e)-FPTAS. 


3.5 Truthful Mechanisms 

This section follows the terminology of | [^ . We define truthful (aka. incentive-compatible) approximation 
mechanisms for our problem. We denote by A C the set of feasible allocations in our problem (NsmCKP 
or MULTl-mDKP). 

Definition 3.2 (Mechanisms). Let V = Vi x • • • x Vn, where Vk is the set of all possible valuations of agent k. 
A mechanism (A, V) is defined by an allocation rule A : V —)• A and a payment rule "P : V —M”. We assume 
that the utility of player k, under the mechanism, when it receives the vector of bids v = {vi, ..., Vn) G V, is 
defined as Uk{v) = Vk{dk{v)) - Pk{v), where A{v) = (di(u),.. .,dn{v)), and V{v) = (pi(u),... ,Pn{v)) 
and Vk denotes the true valuation of player k. 

Namely, a mechanism defines an allocation rule and payment scheme, and the utility of a player is defined 
as the difference between her valuation over her allocated demand and her payment. 
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Definition 3.3 (Truthful Mechanisms). A mechanism is said to be truthful if for all k and all £ Vfc, and 
v-k G V-k, it guarantees that Uk{vk,v-k) > Uk{vk,v-k). 

Namely, the utility of any player is maximized, when she reports the true valuation. 

Definition 3.4 (Social Efficiency). A mechanism is said to be a-socially efficient if for any G V, it returns an 
allocation d G such that the total valuation (also called social welfare) obtained is at least an a-fraction 
of the optimum: r^(d) > a ■ OPT. 

As in ||^[T^[^, our truthful mechanisms are based on using WCG payments with Maximal-in-Range 
(MIR) allocation rules: 

Definition 3.5 (MIR). An allocation rule A. : V ^ X is an MIR, if there is a range TZ (f X, such that for any 
u G V, Al(p) G argmaXjjgT^ u(d). 

Namely, A is an MIR if it maximizes the social welfare over a fixed (declaration-independent) range TZ of 
feasible allocations. If is well-known (and also easy fo prove by a VCG-based argumenf) fhaf an MIR, com¬ 
bined wifh VCG paymenfs (compufed with respect to range TZ), yields a truthful mechanism. If, additionally, 
the range TZ satisfies: maxde 7 ^ p(d) > a • maxdgA" p(d), then such a mechanism is also a-socially efficient. 

Finally a mechanism is computationally efficient if it can be implemented in polynomial time (in the size 
of the input). 

3.6 Application to Power Allocation in Alternating Current Electric Systems 

Conventionally, the demands in (AC) alternating current electric systems are represented by active power in 
positive real numbers and reactive power in (positive or negative) imaginary real numbers, which are complex 
numbers in the first and fourth quadrants of the complex plane. We note that our problem is invariant, when 
the arguments of all demands are shifted by the same angle. For convenience of notation, we assume the 
demands are rotated by 90 degrees unless all demands are entirely in the first quadrant of the complex plane|^ 
See Fig. for a pictorial illustration. Notice that the axis labels are swapped after the rotation, the real axis 
indicates the reactive power while the imaginary indicates the active power. 



Figure 2: Each vector represents a power demand d. The figure shows thaf all demands are rofafed by 90°. 

CKP is a simplified model of fhe real-world AC elecfric systems, which considers fhe capacity consfrainf 
of a single link (e.g., a single boffleneck). There are pracfical scenarios, where fhe single-link capacify con- 
sfrainf is critical. For example, in a microgrid, there is usually a single transmission line connecting the main 

^ Note that it is customary in the power engineering literature to assume that the demands lie in the first and fourth quadrants of 
the complex plane. However, for convenience of presentation, we prefer to work in the first and second quadrants. For instance, if 
all the demands are capacitive (i.e., lie in the fourth quadrant), then rotation by 90 degrees allows us to assume that the all numbers 
involved are non-negative, a property which is necessary for obtaining a PTAS. 
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grid and the microgrid. In such a setting, our model can capture power allocation considering the capacity of 
transmission bottleneck. We also remark that our single-link model is fundamental to general network setting 
with multiple capacitated links. A thorough understanding of the single-link case can pave the way to solv¬ 
ing the multi-link case. As a follow-up study, our recent paper | [T5| considers non-strategic power allocation 
of power flow of inelastic demands in a power network, which is based on some of the fundamental results 
obtained in this paper. 

In the conventional models of AC electric systems, there are multiple operating constraints, in addition to 
capacity constraint. One example is the nodal voltage constraint at customers. In our recent study of event- 
based demand response management in microgrids | |T^ , we evaluate the changes of nodal voltage in response 
to power allocation decisions of customers. We observe that the voltage is less sensitive to the decisions of 
power allocation. Hence, CKP is a suitable model for approximating the power allocation in such a setting. 

We remark that the simplified DistFlow model |[^[3 201 is a well-known approximation model of power 
flows by ignoring the loss terms from the formulation. In fact, CKP is equivalent to the simplified DistFlow 
model on a capacitated single link topology. 

Naturally, MULTiCKP models combinatorial power auctions in AC electric systems. Each customer k 
declares to the utility company a set of demands that represents her preferences among different alter¬ 
natives of load profiles, and a valuation function Vk over Dk- The valuation Vk{(l) represents the amount 
customer k is welling to pay if her load profile d G is satisfied. If customer k wants to bid for a load 
that represents multiple appliances at once, she can include her corresponding vector sum to the set Dk as an 
additional preference. The monotonicity of Vk with respect to the partial order “A” implies the free disposal 
of extra supplied power (see Fig. [T]). Larger active power (imaginary component) should have at least the 
value of smaller active power. Conventionally, capacitive loads have negative reactive power, while inductive 
loads have positive reactive power. Monotonicity implies that extra supplied reactive power of the same type 
(capacitive or inductive) should have at least the same valuation. Indeed, customers can define constant valu¬ 
ation for different values of reactive power. We remark that all demands in the set Dk are assumed to be either 
in the first quadrant or the second, but not in both. This in fact implies the load profiles of each customer 
are either only capacitive, or only inductive. Such mild assumption is actually needed in order to apply the 
(1,1-h e)-FPTAS in Sec.|^ 


4 Hardness of Power Allocation in AC Electric Systems 

In this section, we present our main hardness result for CKP, which depends on the maximum angle 4> the 
demands make with the positive real axis. When (j) £ [| -|- 5, vr], we show that the problem is inapproximable 
within any polynomial factor if we do not allow a violation of Constraint ([^. Moreover, when (j) approaches 
TT, there is no (a, /?)-approximation, for any a and /3 with polynomial bit length. Our hardness results indicate 
that the approximability of the problem CKP differs depending on maximum argument of any demand cj). 
This insight suggests to study different techniques in the later sections to achieve the best approximation 
result possible for each case. 


Theorem 4.1, Unless P=NP, for any (5 > 0 and 6' > 0 

(i) there is no (a, l)-approximation for CKP[^ -|- 6, vr] where a, 5 have polynomial length. 

(ii) there is no (a, 13 ^approximation for CKP[7r — 6', vr], where a and f3 have polynomial length, and 5' is 
exponentially small in n. 


Remark. In fact, these hardness results hold even if we assume that all demands are on the real line, 
except one demand dm+i such that arg(dm_|_i) = | -|- 0, for some 6 £ [<), |] (see Fig.|^. Note that the trivial 




approximation of picking the user with the highest feasible value does not obtain - approximation, because 
we allow demands to have both positive and negative real parts, which can cancel each other. 



Proof. We present a reduction from the (weakly) NP-hard subset sum problem (SubSum): given an instance 
I, a set of positive integers A = {oi,..., am} and a positive integer B, does there exist a subset of A that 
sums-up to exactly B? We assume that B is not polynomial in m, otherwise the problem can be easily solved 
in polynomial time by dynamic programming. 

We construct an instance I' of CKP[| +0, | + 0] for each instance I of SubSum such that if SubSum(/) 
is a “yes” instance then the optimum value of CKP[| + 0, | + 6], denoted by Opt, is at least 1; and if 
SubSum(/) is a “no” instance, then Opt < a even when Cons, (j^ can be violated by /3. 

We define n = m + 1 demands: for each a^, k = 1 ,..., m, define a demand = a^, and an additional 
demand 

dm+i — —B + iB cot 9. 

For all fc = 1,..., m, let valuation Vk = and Vm+i — 1- We let 

C = Bcote. 


We prove the first direction, assuming SubSum(/) is feasible. Namely, = -S, where x G 

{0,1}™ is a solution vector of SubSum. Construct a solution x G {0,1}™+^ of CKP such that 


Xk = 


Xk 

1 


if A; = 1,..., m 
A k = m + 1. 


In fact, this is a feasible solution that satisfies Constraint (j^: using — i? = 0, we get 

/ ra \ ‘^ / m ^^2 

( ^ d^Xk + dm+1 I + ( ^ d\xk + d} 


^m+1 


\k=l 


\k=l 


= d^Xk — B j + B^ cot^ 6 
= B^ cot‘^e = c^. 


Since Vm+i = 1, the total value of such solution v{x) > 1, which implies that Opt is at least 1. 

Conversely, assume that Opt > a. Let x* G {0, l}™+i be an optimal solution that may violate Cons. ([^ 
by I3. Since user m + 1 has valuation Vm+i = while the rest of users valuations total to less than a: 
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< OL, user m + 1 must be included in the optimum. Therefore, substituting in Constraint ( [T^ , 

/ m \ ^ 

dfxl -b\ +B‘^ cot^ e < 


\k=l 


gives 


'Y — Bj < j3‘^C‘^ — B^ cot^ 6 

= B^ cot2 0(^2 _ 


^fc=i 


By the integrality of the Oj’s, 


Yakxl = B ^ \ Yakxl-B\ < 1 


( 12 ) 


(13) 


k=l 


k=l 


In other words, SubSum is feasible if and only if the absolute difference | YlT=i (kkxl. — B\ < 1. This implies, 
SubSum(/) is feasible when the R.H.S. of Eqn. (12i is strictly less than 1. When (3 = 1, R.H.S. of Eqn. ( [T^ 
is zero, and we complete the second direction and hence, the proof of part (i) of the theorem. 

Eor large enough 6, the R.H.S. of Eqn. (12 1 is strictly less than 1: 

B^ cot^ e{l3‘^ - 1) < 1 


This implies, 6 > tan ^ . By Eqn. (13), SubSum is feasible which completes the second 

direction and establishes part (ii) of the theorem. 

□ 


5 A Truthful PTAS for MULTiCKP[0, f - (5] 

In this section, we present our truthful PTAS for MULTlCKP[0, f — 5]. This PTAS invokes a truthful PTAS 
for MULTl-mDKP as a subroutine. Problem MULTl-mDKP was shown in | [T^ to have a (1 — e)-socially 
efficient truthful PTAS in the setting of multi-unit auctions with a few distinct goods, based on generalizing the 
result for the case m = 1 in We explain this result hrst in our setting, and then use it in Sections |S!2] and 


our PTAS works even for = 0. However, we are only able to make it truthful for any given, but arbitrarily 
small, constant 5 > 0. Removing this technical assumption is an interesting open question. 


5.3 


to derive a truthful PTAS for MULTlCKP[0, 5 — <5]. We remark that, without the truthfulness requirement. 


5.1 A Truthful PTAS for MULTl-mDKP 


We present a truthful PTAS for MULTl-mDKP that will be needed in Sec. |5.2p.3| below. This result is a 
slight generalization of | [T8| to accommodate real-valued demand vectors instead of integer-valued. 

Eet c = (c^,..., c™) be the capacity vector, and for any d G 2? C M™, write dk = {d \,..., d^). Eor any 
subset of users N C J\f and a partial selection of demands d = (d^ £ T> : k £ N), such that Ylk^N ^k < c, 
dehne the vector bj^ ^ = (ftjy d> • • • > d) ^ follows 


T i _ ^ 0^ Y^keN ^k 
^N,d (n-|A|)2 • 


(14) 
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Following 18 2^, we consider a restricted range of allocations defined as follows: 


5 = 


U 


^Ar,d) 


( 15 ) 


AfCA/', d = (dj.: keJV): |Af|<^, 

d^ev VkeN 


where, for a set Af C AA and a partial selection of demands d = G 2? : fc G N), 


c _ A 
‘^N,d — 


{(di,...,d„)Gl?"| 5] dk < c,dk = dk Vfc G 




H 3rl G Z+ s.t. 4 = 4 • 4 a and Y,rl<{n- |iV|)2 }. 

k^N 

Note that the range 5 does not depend on the declarations Di,, Dn- The following two lemmas establish 
that the range 5 is a good approximation of the set of all feasible allocations and that it can be optimized over 
in polynomial time. The first lemma is essentially a generalization of a similar one for multi-unit auctions 
in Q, with the simplifying difference that we do not insist here on demands to be integral. The second lemma 
is also a generalization of a similar result in |[8j, which was stated for the multi-unit auctions with a few 
distinct goods in 1181. For completeness, we give the proofs in the appendix. 

Lemma 5.1 ( |[^). max^g^ n(d) > (1 — e)OPT. 


Lemma 5.2 (|8 18|). We canfindd* G argmax^jg^ n(d) using dynamic programming in time 


0(m/e) 


It follows that an allocation rule defined as an MIR over range S yields a (1 — e)-socially efficienl frulhful 
mechanism for MULTl-mDKP. 


5.2 A PTAS for MULTlCKP[0, f ] 

We now apply fhe resulf in fhe previous section fo fhe mulli-minded complex-demand knapsack problem, 
when all agenfs are restricfed fo reporf fheir demands in fhe posifive quadranf. We begin firsf by presenting a 
PTAS wifhouf sfrafegic considerations; fhen if is shown in fhe nexf section how fo use fhis PTAS wifhin fhe 
aforemenfioned framework of MIR’s fo obfain a frufhful mechanism. 

Overview of Technique. As we will see in Sec. if is possible fo obfain a (1,1 -|- e)-approximafion by 
a reduction fo fhe MULTI-2DKP problem. To gel a heifer resulf wifhouf violaling fhe constrainl, we reduce 
MultiCKP[ 0, |] inslance fo MULTl-mDKP. We note lhal allhough Ihere is a PTAS for MULTl-mDKP wilh 
conslanl m 0, such a PTAS cannol be direclly applied fo MULTlCKP[0, |] by polygonizing fhe circular 
feasible region for MULTlCKP[0, |], because one can show fhal such an approximation ratio is al leasl a 
conslanl faclor. This is fhe case, for inslance, if fhe optimal solufion consisfs of a few large (in magnilude) 
demands logelher wilh many small demands, and if is nol clear al whal level of accuracy we should polygonize 
fhe region fo be able fo caplure Ihese small demands. To overcome Ibis difficully, we have to first guess the 
large demands, then we construct a grid (or a lattice) on the remaining part of the circular region, defining a 
polygonal region in which we fry fo pack fhe maximum-ulilify sel of demands. The lafler problem is easily 
seen fo be a special case of fhe MULTl-mDKP problem. The main challenge is fo choose fhe granularity of 
fhe grid small enough fo well-approximale fhe oplimal, bul also large enough so lhal fhe number of sides of 
the polygon, and hence m is a constant only depending on 1 /e. 

In this section we assume that aig{d) < |, that is, > 0 and 4 > 0 for all d £ T>. Without loss of 
generality, we assume e <i where ^ G Z_|_. For an integer i G Z+, let Ci{i) and € 2 ( 1 ), respectively, denote 
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the sets of all vertical and all horizontal lines in the complex plane that are at (non-negative) distances, from 
the real and imaginary axes, which are integer multiples of S, that is. 


- {a; + iy G c I X = ^, A G Z+,y e M}, 
AC 

J- 2 U) = {x + iy GC\y = ^, XgZ+,x gR}. 


Given a feasible set of vectors T C 2? to MULTlCKP[0, |] (that is, | Yld&T define dr = Yld&T 

and let 

tur — \J— Re(dT)^ — Ini(rfr)) Wt — y/C^ — Im(dr)^ — Re((ir)- ( 16 ) 

Namely, (resp. is the horizontal (resp. vertical) distance between (It and the boundary of the con¬ 
straint disk (see Fig.j^for an illustration). 

Let pi{T) and P 2 {T) be the smallest integers such that 

C/2Pi(^)<^andC/2^^(^)<^. 

The set of lines in Ci{pi{T)) U C 2 {p 2 {T)) define a grid on the feasible region at “vertical and horizontal 
levels” pi{T) and P 2 {T), respectively. We observe that if we increase d^, then decreases as well as the 
real granularity of the grid, and hence pi (T) becomes larger. Similar observation holds when we increase d^ 
Let Ai(r) and \ 2 {T) be the largest integers such that 


dr > 


Ai(r)c 


and dr > 


\2{T)C 

2P2(T') ’ 


and ^r G C be the intersection of the two lines corresponding to Ai(r) and A 2 (T): 


ZT = {x + iy gC\x = } n{x + iyGC\y= 

Given zt, we define four points in the complex plane (V’^r; V'r; V’^r) that 

= ( 0 , v'C2-Re(zT)2), V’r = (Re(zr), 
tp'r = (\/G2 - Im(zT)^,o), 'i/’r = (\/G2 - Im(zr)^, Im(2;T)) ■ 

Let TZt be the part of the feasible region dominating zt- 


TZt = {x -|- iy G C : |x -|- iy| < C, x > Re(2;r), y > Im(zr)}, 

and Rr(e) be the set of intersection point^between the grid lines in £i(/3i(T))u£2(p2(R)) and the boundary 
of T^r’ 

Pr(e) = {^G 7 ^r : |z| = C} n (£i(pi(r)) U £2(P2(T))). 

The convex hull of the set of points Rr(e) U {V^'r) V’r; V'r; defines a polygonized region, which we 

denote by 'Pr(e) and its size (number of sides) by mr(e) (see Fig.|^for an illustration). 

Lemma 5.3. mr(e) < ^ + 3. 

^For simplicity of presentation, we will ignore the issue of finite precision needed to represent intermediate calculations (such as 
the square roots above, or the intersection points of the lines of the grid with the boundary of the circle). 
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Proof. Let x be the horizontal distance from zt to the boundary of the circle (with center 0 and radius C). 
Then, by the definition of pi{T), which implies that ^ On the other hand, by 

the definition of zt. 

It follows from the above inequalities that the number of vertical grid lines (at level pi(r)) between zt and 
the boundary of the circle is at most 


C/2/’i(^) 


+ 1 < 


2Pim (l + I) y;. 


c 


+ 1 < 


8{l + i) 


+ 1 < 


9 

e 


Similarly, we can show that the number of horizontal grid lines between zt and the boundary of the circle is 
at most Adding the three other points {'f'j'., 0} gives the claim. □ 



Figure 4: We illustrate the region TZt by the shaded area and Pt{^) by the black dots on the arc of the circle. 
The zoomed area highlights that (resp. wlp) is defined by dx not zt- 

Definition 5.4. Consider a subset of users N C M and a feasible set T = {dk '■ k G N} to MULTlCKP[0, |]. 
We define an approximate problem (PCZ^j by polygonizing MULTlCKP[0, |].' 

(PGZt) max ^ 0^(4) 

k£jV 

s.t. 'Y^dke Vt{() 

keJV 

dk = dk, yk € N 
dk £V, yk£ M\N. 

Given fwo complex numbers p. and o, we denote fhe projecfion of /r on o by Pjj,(/r) = j^{p^u^ + p}i'^). 

Given fhe convex hull VxiPj, we define a sef of mT{e) vectors {erf}, each of which is perpendicular to each 
boundary edge of Vxie) and sfarfing af fhe origin (see Fig.j^for an illusfralion). 

Definition 5.5. Consider a subset of users N C M and a feasible set T = {dk : k G N} to MULTlCKP[0, |]. 
We define a MULTl-mDKP problem based on {cr^}.' 
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Im 



Figure 5: Each cj\, is a vector (starting at the origin) perpendicular to each boundary edge of Vt{^)- 


(MULTl-mDKP{(T^}) max Vk{dk) (17) 

k£jV 

|Pja^f’(4)| < krl, Vi = 1,... ,mr(e) (18) 

kej\f 

dk = dk, yk € N (19) 

4 e 77, V/c G M\N. (20) 


One can see that MULTl-mDKP{(Ty} is an instance o/MULTl-mDKP (defined in Sec. 3.3) over the set of 
users N\N, m = mT{e), d\ = |Pj,,T(4)|, and c* = \afi\ - YlkeN\^iaf(dk)\. 

Lemma 5.6. Given a feasible set T to MULTlCKP[0, |], PGZ^ and MULTl-mDKP{cjy} are equivalent. 


Lemma 5.6 follows straightforwardly from the convexity of the polygon Vrie)- 

Our PTAS for MUTLlCKP[0, |] is described in Algorithm MultiCKP-PTAS, which enumerates every 
subset partial selection T of at most 7 demands, then finds a near optimal allocation for each polygonized 
region Vt{() using the PTAS of MULTl-mDKP from Section 5.1 which we denote by MULTl-mDKP- 
PTAS[]. 


Algorithm 1 MultiCKP-PTAS ({ufc, Dk}k&^f, C, e) 

Require: Users’ multi-minded valuations {vk, Dk}ke^f', capacity C; accuracy parameter e 
Ensure: (1 — 3e)-allocation [di, ..., 4) to MULTlCKP[0, |] 

1 : ( 4 ,..., 4 ) ( 0 , ..., 0 ) 

2: for each subset N C Af and each subset T = {4 £ Dk ■ k G N} of size at most 7 s.t. | YldeT d\ < C 

do 

3: Set dx ^ Yld&T define fhe corresponding vectors {erf} 

4: Obfain (di,..., 4) MULTl-mDKP-PTAS [MULTl-mDKP{iT^}] wifhin accuracy e 

5: if Y.kf{^k{dk)< Efc '^k{dk) then 

6: • • • ) dn) i {di ,. . . , dji) 

7: end if 

8: end for 

9: return (di,..., 4) 
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Theorems.?. For any e > 0, Algorithm MULTlCKP-PTASyz«<i5a (1—3e, 1)-approximation to MULTlCKP[0, ^]. 
The running time of the algorithm is ||J^ 



Vk{dk) > (1 - e) Vk{dk) > (1 - 3e)OPT. 

k k 

We give an explicit construction of the allocation (di,..., (i„) in Algorithmic thus completing the proof by 
Lemma 15^ □ 


Lemma 5.8. Consider a feasible allocation d = (di ,..., d„) to MULTlCKP[0, ^]. Then we can find a set 
T C {di,..., dn} and construct an allocation d = (di,... , dn), such that |r| < ^ and d is a feasible 
solution to VGXt and r'(d) > (1 — 2e)n(d). 


Proof In Algorithmic let I and be the values of ^ and Ti at the end of the repeat-until loop (line|Cl. 

The basic idea of Algorithm |C is that we first construct a nested sequence of sets of demands To C 
Ti <Z ... <Z Ti, such that a demand is included in each iteration if it has either a large real component or a 
large imaginary component. The iteration proceeds until a sufficiently large number of demands have been 
summed up (namely, \Ti\ > ^), or no demands with large components remain. At the end of the iteration, if 
the condition in line[^holds, then S = T^, i.e., the whole set S can be packed within the polygonized region 
'PT^(e). Otherwise, we find a subsef of S fhaf is feasible fo PGZ^y. 

To do so, we parfifion S\Ti info groups (possibly one group), each having a large componenf along eifher 
fhe real or fhe imaginary axes, wifh respecf fo fhe boundaries of fhe region Then removing fhe group 
wifh smallesf value among fhese, or removing one of fhe large demands wifh smallesf valuation will ensure 
fhaf remaining demands have a large value and can be packed within ^^^(e). 

We then have to consider two cases (line|^: (i) \Ti\ becomes at least ^ or (ii) S'|*' USj = 0. For case (i). 
We reduce the size of Tg by taking the first - demands, and then we combine the remaining demands in S\Ti 
into a group Vi. We show next that removing any one demand G will make S\{dk} a feasible solution 
to PGZ'p-yjrf^}. Since dx^ F dx^, for i > C , the lengths and are monotone decreasing for ^ = 1,2,... 

Indeed, w.l.o.g., let I be such that d^ G (also d^ 0 Tf) and let T = Tf\{dk}. Then, implies 


that 


C C 


< 


ew. 


< d^ 


k ■ 


Hence, d^ is larger than the width of the grid’s cell Therefore, is a feasible solution to PGZ-r 

because the removed demand exceeds the width of the grid (see Fig.j^for an illustration). 


For case (ii), we can apply Lemma 5.9 below to partition 5'\T^ into at least 7 — 1 groups {Vi,..., I 4 }, 


where each group Vj has a large total component along either the real or the imaginary axes (precisely, greater 
than or ^wJp_ respectively). This implies that removing any group Vj will make T^C{Jj/^j Vj' a feasible 
solution to PGZx^. 
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Figure 6: is larger then the grid’s width. If dk is removed, the remaining set of demands is feasible to 

PGZt- 



Figure 7: Demands {dk \ T G AA\T} are partitioned into {Vi}. 


To conclude, there are either (i) at least ^ demands in |T^|, or (ii) Sf U Sj = 0. We define S' by deleting 


a minimum valuation demand or group of demands from S (lines 


23 


and 


261. Then, we set dk = dk if k ^ 


and dk = 0 if k ^ S'. Hence, in case (i), u(d) > (1 — e)OPT, and in case (ii), u(d) > (1 — ^)Opt 
■ Opt > (1 - 2e)OPT. 


Lemma 5.9. Consider a set of demands S CT) and TVS, such that 

1. S is feasible solution to MULTlCKP[0, ^], but S is not a feasible solution to POZy 

2. dP^ < and d^ < for all d G S\T. 

Then there exists a partition {Vi,... ,Vh} of S\T such that 

1. either (i) Yhd&Vj all j = I,..., h, 

2. or(ii) ^ {wlpfor allj = l,..., h. 

where /i G [2 — 1, |). 

Proof. First, we define k = 'f2d&s\T If 5' is a feasible solufion fo MULTlCKP[0, |], buf S is nof a feasible 
solufion fo PGZt, then at least one of the following two conditions must hold: either (i) or (ii) 

p > (see Fig M for an illustration). 
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Algorithm 2 CONSTRUCT({T;fc}fcg;V '5 d, C, e) 

Require: Users’ valuations {t'fc}fceA/'i ^ feasible allocation d = {di,dn)', capacity C; 
accuracy parameter e 

Ensure: A set of demands T C {di,..., and a feasible allocation d 
Initialize: S •«— {di,..., d„}; d = d; ^ ^ 0; ■(— 0; ^ 0 

t> Find a subset of large demands 

1: repeat 

2: i^i+l 

3: dr^ ^ Y,d£Te-i ^ 

4: ^/C^^^^Re(dT^ - Im(dTj; ^ - Im(dTj^ - Re(dTj 

5: ^ {d E 5 \ T, I d^^ > IS] ^ {d E S \ d^ > 

6: ^ U Sf U S\ 

7: r]i_i + ^ 

8: until IT^I > f or U S'] = 0 or S\Ti = 0 
9: Yld£S\Te ^ 

10: if S\Ti = 0 or r]£ + K G Vt^ (e) then 

11: return (T^jd) 

12 : else 

0 Find a subset S' G S that is feasible to PCZ-p^ 

13: if IT^I > f then 

14: Ti G- the set of the first - elements added to Ti 

15: h^l\Vi^ S\T£ 

16: else 

17: Find a partition Ui,..., 14 over S\T£ such that either 

(i) Edeu for all j = 1,, h, or 

(ii) Edeu for all j = 1,..., h, 

where h is defined in Lemma l5^ 

18: end if 

19: Pick k E argmin{nfc(dfc) | d^ E T^} 

20: Pick j E argmin{X;A::d;,eU \j = ^,---,h} 

21: if v^{d^) < Ek:d,&v. Mdk) then 

22: ^ 0 

23: return (T£\{d^},d) 

24: 

25: dfc •(— 0 for all A: : dfc E IF 

26: return {Ti, d) 

27: end if 

28: end if 


Without loss of generality, we assume case (i). Let us pack consecutive demands d E S\T into batches 
such that the sum in each batch has a real component of length in ^w^]. More precisely, we fix an 
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order on S\T = {di, ..., dr}, and find indices 1 = ki < k 2 < ■ ■ ■ < kh' < kh'+i = r + 1 such that 


kt+i-l 

fort = 

k=k£ 
ki+i 

and ctfc > for t = 1,..., — 1. 

k=ki 


( 21 ) 

( 22 ) 


It follows from Eqn. (22 1 that Ylk=ke ^ > t^T for t = 1,..., /i' — 1, since < f It also follows 


that - < h' < - + 1, since summing Eqn. (21 1 for £ = 1,..., /i' yields 


h' —1 r 

h'w^ > ^ Yj df = Y^k = > 

i=l k=kj( k=l 




R 


(23) 


The last inequality follows from our assumption. 


Similarly, summing Eqn. (22 1 for £ = 1,..., /i' — 1 yields 


h'—l 


{h' - l)^w^ <Y Y 


£=1 k=ki 


k=l 


where the last inequality is derived by the feasibility of S. Setting Vf = ..., dfc^^^_i}, for i = 

1,..., /i' — 2, = {dk^,_^ ■,■■■, dr}, and h = h' — 1 satisfies the claim of the Eemma. □ 

5.3 Making the PTAS Truthful 


In this section, we make the PTAS, presented in Sec. 5.2 truthful. One technical difficulty that arises in this 
case is that the polygons 'P'r(e) defined by a guessed initial sets T are not monotone w.r.t. the set of demands 
in T, that is, if we obtain T' from T by increasing one of the demands from dk to d'^ dk, then it could 
be the case that Vt{^) 2 'Pt'{^)- Hence, it is possible to manipulate the algorithm by a selfish user in T 
who untruthfully increases his demand to change his allocation to become a winner. To handle this issue, 
we will show that the number of possible polygons that arise from such a selfish user, misreporting his true 
demand set, and can possibly change the outcome, is only a constant in e and <5 (recall that we consider 
MultiCKP[ 0, I — d]). Thus, it would be enough to consider only all such polygons arising from the reported 
demand set. Eigure [^provides a pictorial illustration. 

Since we assume that arg(d) G [0, | — d], for all d G |J^ D^, we may assume further by performing a 
rotation that any such vector d satisfies arg(d) G [f, f — |]- For convenience, we continue to denote the new 
demand sets by D^, and redefine the valuation functions in terms of these rotated sets. By this assumption. 


d dl / d^ 

tan - < —5- < tan - 

2 - “ V 2 


for any TCP. 


(24) 


We may also assume, by scaling e by 2/(1 + 2 cot^ |) if necessary, that 


e < 


1 + 2 cot^ i 


( 25 ) 


Now we state an important lemma that will be used to prove our main result. 
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Figure 8: We consider two polygonized regions TZt and TZt', one defined by dx and the other by dx' (which is 
obtained by slightly increasing one demand in T). The hgure illustrates the case that such a slight increase can 
possibly change the granularity of the grid such that an infeasible solution in Vt{() would become feasible in 


Lemma 5,10. Let T,T' CP be such that dx ^ dr'- Consider a vector k G C such that dx' + n € Vx'{C). 
Then either (i) dx + n ^ ‘Pr(e). or(ii) pi{T') < pi{T) + 1 and P 2 {T') < P 2 {T) + 1- 

Proof. Suppose that dx + n ^ Vx{C). Since dx' ^ dx, it also holds that dx' + k 0 Vx{C). This implies 

that both dx + K and dx' + rt he within the same grid cell at vertical and horizontal levels pi{T) and P 2 {T), 

respectively, because otherwise \dx + K + > C and dr + « + 0 ^r'(e) 


which is a contradiction. Hence d^, — d^ < ^ 

From the definition (16 1 of w^, we have 

(w^ + d^f = C^ - {d^xf 


< and dip, — dL < 


c 


T - 2P2(t) - 


and {w^, + df;/)^ = — (dlp,)'^ 

Subtracting the above equations from each other obtains 

{wx + dy)^ — {w^, + dy/)^ = (dy/)^ — (dy)^ 

^^(^Wx T d^ — w^/ — d^/)(?n^ -t- d^ T w^, d^/) = (d^/ — d^)(d^/ T d^) 


^^Wx “b dx — Wrp, — drp, — (^drpf — d-j-) 




„R 


R 


d^/ “b d}p 


w, 


R 


“b dip “b W^f “b dip, 


^X' 


Therefore, 


Wx 


— Wxf “b d^/ — d^ “b iwx “b d^ — wl^, — d^/) 
= wl^, + d^, — d^ + [dip, — dip) ^ 


dip, “b dip 


d^, + w^, + dp + w. 


S Wp, H - 


R 


+ 


ewt, 


rp f drp, ~b drp 


4 \ dip, + dl^ 


R emS 
< wf, + —^ 


Wtp 

1 H-p • 

wr 


1 


p tan 2 


(26) 


(27) 


(28) 


where we use Eqn. (261 in Eqn. (27), and use Eqn. (24i in the last inequality. We can upper-bound Wp/wj 
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by 2 / tan | also using (241 as follows: 











< 


1 + 


~C 





d5 


) + 7 ?^ tan 7 


C 


The latter quantity is bounded by /(I) = ^ r , since the function /(a) = 


t <5’ 

tan 2 


\/l—a^+atan f 


IS monotone 


increasing in a G [0,1]. Using this bound in ([27|) and rearranging terms, we get 


{ 1 - |(1 + 2cot2 


(29) 


by our assumption (251 on e. From (291 and 


T' 


< 


pi(T) + 1. Similarly, we have P 2 {T') < P 2 {T) + 1. 
We now state our main result for this section. 


c 

2Pi(t') 


and 


c 

2PpT) 


< follows that pi{T') < 

□ 


Theorem5.il. For any e,5 > Othereisa {1—3e)-socially efficient truthful mechanism for MvuiCKP[0, ^ — 
(5]. The running time is ||J^ v ' 

Proof It suffices to define a declaration-independenf range S of feasible allocations, such fhaf maxjg^ n(d) > 
(1 — 3e) • Opt, and we can optimize over S in fhe sfafed time. 

For T C X), lef G(T) be fhe sef of vectors in C defined by fhe union of and 

(a) fhe (componenf-wise) minimal grid poinfs z G TZt, such fhaf z = iiFi £2 for some ii G £1 (pi (T) + 1) 
and X 2 G C, 2 {p 2 {T) + 1), and eifher pi{{z}) = pi{T) + 1 or P 2 {{z}) = P 2 {T) + 1, buf nof bofh; and 

(b) fhe (componenf-wise) minimal grid poinfs 2; G TZt, such fhaf z = iiCi £2 for some G G £1 {pi (T) -|- 1) 
and £2 G C 2 {p 2 {T) + 1), and pi{{z]) = pi(T) -|- 1 and P 2 {{z}) = P 2 {T) + 1- 

Nofefhaf |G(r)| = 0(f). Wemenfionfhafifispossiblefoenumerafeoverall£i(pi(r)-|-l)n£ 2 (/ 02 (T')-|- 
1) poinfs (i.e., fhe infersecfion poinfs of fhe doffed grid lines in Figure]^, insfead we choose to only enumerafe 
over a smaller subsef of Parefo minimal poinfs fhaf are defined by (a) and (b). Figure gives a pictorial 
example of fhese minimal poinfs. 

For convenience of nofafion, lef us fix fwo subsefs Pi, X ?2 U V. For 2 G G{T), lef us denote by Sz(V 2 ) 
fhe range of feasible allocafions defined as in ( [T5] ) wifh respecf to fhe MULTl-mDKP problem wifh consfrainfs 
(fT8])-([20l), when 

(I) T is replaced by T U {z — dr} (and hence, z is used to define fhe polygon 'PT(e)); 

(II) we add an additional “dummy” user n -|- 1 to A( with valuation Vn+i{d) = 0 for all d G P, such that 
the vector z — dr allocated to this user; and 

(III) the set of vectors in Af\N is chosen from P 2 . 
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Figure 9: The figure illustrates the points of the set G{T). The black solid lines are Ci{pi{T)) and C 2 {p 2 {T)) 
on the vertical and horizontal directions respectively; the gray dotted lines (with some overlapping with the 
black lines) correspond to Ci{pi{T) + 1) and C 2 {p 2 {T) + 1)- The red square points correspond to condition 
(a), and the blue diamond points correspond to condition (b). 


Then we define the range S{'Di,T> 2 ) as the union: 


5(2?i,P2)= U 

TCDi: |T|<i 


.zeG(T) 


By Lemmas 


5.1 


and 


that we can efficient 


5.8 


we have maXdg 5 ( 25 ^D) u(d) > (1 — 3e)OPT (since dT G G{T)). It remains to argue 


y optimize over S{'D, V). Using Lemma 5.10[ we argue that we can solve the optimiza¬ 


tion problem over S{V, V) assuming that V = |J^ that is, maXdg 5 (x)^x)) u(d) = max^g^dj^ 


^fciUfe ^fc) 


v{d). 


One direction “>” is obvious; so let us show that maxdg 5 (x)^x)) u(d) < max^g^dj^ D^) ^(d). 

Suppose that d* = (dj,..., d*) is an optimal allocation over S{V, T>), but such that d* G for some 
z' G G{T"), T" C V, and T" % |J^ Dk- Then let us show that there is a set T C |J^ D^, z G G{T), and 
d G Sz(T>), such that u(d) = u(d*). 

Define an allocation d as follows: Let N = {k : d^ G T"}; for each /c G A^, we choose dk G Dk such 
that dk ^ d|, and Vk{dk) = Vk{d\), and we keep dfc = d|, if fc 0 N. Let us apply the statement of the lemma 
withT = {c 4 : A: G N],T' = r"U{z'-dT"}, and k = Yuk-.k^Nyj{n+i} ^k- If (i) holds then dr + k G 
and therefore we have 

max ufd) = max ^(d). (30) 

d&s{v,v) de-S(Ufe^fc,i’) 

On the other hand, if (ii) holds, then pi{T') G {pi(T), pi(r) + 1} and P 2 {T') G {p 2 {T)^ p 2 {T) + 1}. In this 
case, if pi{T') = piiT) and P 2 {T') = P 2 {T) then VT'iG) C Vt{(^) (since d^ ^ dx'), in contradiction that 
(i) does not hold; otherwise, there is a point z G G{T) such that z ^ z', pi{T U {z — d^}) = pi{T') and 


P 2 {T U {z - dr}) = P 2 {T'). Then z + k^z' + kG C T’ru{z-dT}(^)’ 

Finally, we note that 


max u(d) = max ^^(d), 

d£S{\Ji^Dk,V) decS(Ufc£>fe,Ufc-Dfe) 


as follows from (the proof of) Lemma [5?^ 


□ 
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6 A Truthful (1,1 + e)-FPTAS for MULTlCKP[0, tt - ^] 

In this section, we present our truthful (1,1 + e)-approximation for MULTlCKP[0, vr — <5] by a reduction to 


MULTI-2DKP. As in 1181, the basic idea is to round off the set of possible demands to obtain a range, by 


which we can optimize over in polynomial time using dynamic programming (to obtain an MIR). 

Let 9 = max{(/) — f, 0}, where cp = arg(d). We assume that tan9 is bounded by an a-priori 

known polynomial P{n) > 1 in n, that is independent of the customers declarations (valuations and demands), 
because the power factors are often limited by certain regulations in practice. 

Let A/+ = {A: G AA I > 0 Vd G Dk} and N- = {/c G AA | d^ < 0 Vd G Dk} be the subsets of users 
with demand sets in the first and second quadrants respectively (recall that we restrict users’ declarations to 
allow such a partition). 

We can upper bound the total projections for any feasible allocation d = (di,..., dn) of demands as 
follows: 


E 4 < c, E —d^ < Ctand, d^ < C(1 + tand), 

keJ\f k£M-{d) fceA4(d) 


(31) 


where A/+(d) = {/c G AA | d^ > 0} C A/+ and Ad (d) = {A: G AA | d^ < 0} C Ad. Define L = 
and for d £ V, define fhe new rounded demand d as follows: 


I if/l A 


/ 

'<p' 

L 

•L + i 

L 





\ 

cP 

L 

•L + i 

'd}~ 

L 


• L, if d^ > 0, 

• L, ofherwise. 


(32) 


For convenience, we will write d = (di,...,d„). Note fhaf by Ibis definition, A/’+(d) = A/’+(d) and 
Ad(d) =AA_(d). 

Consider an optimal allocation d* = (dj,..., d*) fo MULTiCKP [0, tt — d]. Lef (and ,^_), (and 
(j_) be fhe respecfive real and imaginary absolufe fofal projections of fhe rounded demands in A/+(d*) (and 
Ad(d*)). Then fhe possible values of (^_ are integral mufiples of L in fhe following ranges: 


2L ,..., 

C+,C- |o,d2L,..., 


■C(l + P(n)) 
L 

'C-P{n) 


■ L 


L 


■ L 


C 

T 


L 


We firsf presenf a (1,1 + 4e)-approximation algorifhm (MultiCKP-FPTAS) for MULTlCKP[0, vr — d]; 
fhen we show how fo implemenf if as an MIR mechanism. 

The basic idea of Algorifhm MultiCKP-FPTAS is fo enumerate fhe guessed fofal projecfions on real and 
imaginary axes for A/+(d*) and A/’_(d*) respectively. We fhen solve fwo separate MULTI-2DKP problems 
(one for each quadranf) fo find subsefs of demands fhaf salisfy fhe individual guessed fofal projections. Buf 
since MULTI-2DKP is generally NP-hard, we need firsf fo round fhe demands fo gef a problem fhaf can be 
solved efficienlly by dynamic programming. We nofe fhaf fhe violafion of fhe opfimal solufion fo fhe rounded 
problem w.r.f. fhe original problem is small in e. 

For any allocation d = (di,..., dn) G D”, lef us for brevify write ^+(d) = EkeM+{d)dk’ ^-(d) = 

SfceN_(d) D(d) = Yhk&M d^k- Then by (32 1 and fhe facl fhaf x <t\j\ < x + tfor any x, t such 

thaf f > 0, we have 


max{r(d) — nL, 0} < T(d) < T(d), 


(33) 
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for all r G {r+, r_, r/}. 


Lemma 6.1. For any feasible allocation d = (di,..., dn) to MULTiCKP [0, vr — <5], we have | Yl,k < 
(1 + 2e)C*. 

Proof. Using ( |^ and ( |3T] ), 

f Z] + f Z 

\k&M / \k&N / 

= T|(d) + (d) - 2r+(d)r_(d) + Tf(d) 

< (r+(d) + nL)‘^ + (T_(d) + nL)^ - 2T+(d)T_(d) + (T/(d) + nL)^ 

= (T+(d) - T_(d))^ + Tj{d) + 2nL(T+(d) + r_(d) + r/(d)) + 3n^L^ 

= fz^fc) +fz^fc) +2nL |4| + ^ 

V/ceA/" / \k^M / \k&M k&M / 

< + 4nL(P(n) + 1)C + Sn^L^ = + AeC^ + 2,e^C‘^/{l + P[n)f 

< C‘^{1 + 4e + Se^) < C‘^{1 + 2ef. 


□ 


The next step is to solve the rounded instances exactly. This can be done with essentially the same 


dynamic program used in the proof of Lemma 5.2 modulo a small modification, which we include here for 
completeness. 

For d G 2?, we use d to denote the vector in C with components and cF = dF We define a 

(new) valuation function by: Vk{d) = Vk{d), for k G A/+, and Vk{d) = Vk{d), for k G A2_. Let further 
V^ij^GV: d^ GA+andd^ gB}, and note that \V\ = 

Assume an arbitrary order on A2 = {1,..., n}. We define a 3D table, with each entry U{k,c) being the 
maximum value obtained from a subset of users {1, 2 ,..., 2} C J\f, each choosing a demand from V, such 
that the chosen demands fit exactly within capacity c G V (i.e., satisfy the capacity constraints as an equation 
along each of the axes). The cells of the table are defined according to the following rules: 


U{l,c) = vi{L ■ c); 

U{k, c) = —00 for all c ^V; 

U{k + l,c) = max{i;fc_|_i(L ■ (f) + U {k, c — F)}. 

c'GD 

The corresponding optimal allocation (F+ or F_) can be mapped to the original range of demands V as 
follows: 


X(l, c) = {(1, di)} where di G Di is 


J fi(di) = vi{L ■ c) and di A L • c, if 1 G M+, 
\ vi{di) = vi{L ■ c) and di A L • c, if 1 G J\f-\ 

Z{k + 1, c) = Z{k, c) U {(fc + 1, dfc+i)} where d^+i G D^+i is 


f '(;fc+i(dfc+i) = r;A;+i(L • d) and dfc+i A L • d, ifk + lGJ\f+, 
\ Vk+iidk+i) = Vk+iiL ■ d) and dfc+i < L ■ d, if A; + 1 G N-, 
where d G argmax^,gg{'i;fc+i(L • c') + f7(fc, c — c)}. 
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Algorithm 3 MultiCKP-FPTAS ({ffc, C, e) 

Require: Users’ multi-minded valuations {vk, Dk}k^j^-, capacity C; accuracy parameter e 
Ensure: a (1,1 -|- 4e)-approximation d = {di, ..., dn) to MULTlCKP[0, tt — <5] 

1 : {di,...,dn) ^ ( 0 , ..., 0 ) 

2: P ^ {f E P ; d^ G A+and d^ G B} 

3: for all E A+, G C-i-, C- G B do 
4: if(e+-e-)^ + (C+ + C-)^ < (l + 2e)2c2then 

5: F+ ^ MULTI-2DKP-EXACT({i;fc,Ufe}fegAA+, 

6: F_ ^ MULTI-2DKP-EXACT({vfe,Ufc}fcgAA_, ^,V) 

7 : {d [,..., d'n) <— F^LI F- 

8: if Efc Vk{d'^) > Efc Vk{dk) then 

9 : {di,...,dn)^ {d[,...,d'J 

10: end if 

11 : end if 

12: end for 

13: return (di,..., dn) 


This table can be filled-up by standard dynamic programming; we denote such a program by MULTI-2DKP- 
Exact[-]. 

The following lemma states that the allocation returned by MultiCKP-FPTAS does not violate the 
capacity constraint by more than a factor of 1 -|- 4e. 


Lemma 6.2. Let d be the allocation returned by MultiCKP-FPTAS. Then \ dk\ < (1 + 4e)C'. 


Proof. As in the proof of Femma 6.1 

2 

+ 




(T+(d) -r_(d))Vr|(d) 

r|(d) -P E(d) - 2r+(d)r_(d) -P T|(d). 


If both T+(d) and r_ (d) are less than nL, then the R.H.S. of (34 1 can be bounded by 


T+(d) -P r^(d) -P r|(d) < T+(d) -P E(d) - 2r+(d)T_(d) -P -P r|(d) 

= (r+(d) - T_(d))^ -P Tj (d) -P 2 n^L^. 


(34) 


(35) 


Otherwise, we bound the R.H.S. of ( [^ by 

T+(d) -P r^(d) - 2(r+(d) - nL)(r_(d) - nL) + Tj{d) 

= (T+(d) - r_(d))2 -P r|(d) -P 2nL(r+(d) -P r_(d)) - 2n^L^. (36) 

Since d = is obtained from feasible solutions Fj^ and F_ to 2DKP-ExACT({ufc, Dk}ke:M+} 

and 2DKP-EXACT({i;fc, i^fclfceA/'- > respectively, and C+, (- satisfy the condition in Step|^ 
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it follows from (341-( 361 that 


J \k&M / 


- f X/+ 2nL y~] \d^\ + 

XkSiM J \k&M / k&M 
= (^+ — + (C+ + C-)^ + 2nL(^+ + ^_) + 2n^L^ 

{l + P{n))C + 2 n^ 


< ( (1 + 2e) C + An- i i 

V n{P{n) + 1) 

< ((1 + 2ef + 4e + 2e‘^) < (1 + AefC"^. 




n^{P{n) + 1)2 


□ 


Theorem 6.3. For any e,d > 0, there is a truthful mechanism for MULTlCKP[0, vr — 5], that returns a 
(1,1 + Ae)-approximation. The running time is polynomial in n, cot 6, and 

Proof We define a declaration-independent range S as follows. For G ^+, G A-, C+) C- £ define 
5c+,5_,a,C-={d=(di,...,d„)GP": ^L.4 = ^+, J;L.4 = C+, 


fceA/V 


fceA/V 


L-df = ^., Y. L-dl = C-}. 

fceA/"- feeW- 


Define further 


5 = 


IJ 5j+,€_,c+,C-- 

(5+-?_)2 + (C++C-)2<(1+2€)2C2 


Using Algorithm MultiCKP-FPTAS, we can optimize over S in time polynomial in n and Thus, it 
remains only to argue that the algorithm returns a (1,1 + 4e)-approximation w.r.t. the original range D”. 
To see this, let , d* G P be the demands allocated in an optimum solution to MULTlCKP[0, vr — d], 

and di,... ,dn G P be the demands allocated by MultiCKP-FPTAS. Then by Lemma [64] the truncated 
optimal allocation {d^,..., d*) feasible with r^ect to a capacity of (1 -|- 2e)C, and thus its projections 
will satisfy the condition in Step got Algorithm]^ It follows that u(d) > v id*) > r;(d*) = Opt, where 
the second inequality follows from the way we round demands and the monotonicity of the valuations. 
Finally, the fact that the solution returned byMULTiCKP-FPTAS violates the capacity constraint by a factor 
of at most (1 -|- 4e) follows readily from Lemma 6.2 □ 


7 Conclusion 

In this paper, we provided truthful mechanisms for an important variant of the knapsack problem with 
complex-valued demands. We gave a truthful PTAS when all demand sets of users lie in the positive quadrant 
(which is attained by limited power factors), and a truthful FPTAS with capacity augmentation when some of 
the demand sets can lie in the second quadrant (which captures the general setting of large power factors). Our 
hardness results show that this is essentially best possible assuming P^^NP Our fundamental results underpin 
a wide class of resource allocation problems arising in smart grid. The complete understanding of truthful 
complex-demand problem paves the way to more sophisticated and efficient mechanism design in future smart 
grids. Recently, this work has been extended to consider scheduling problems 112 T3]T7|. 
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APPENDIX 

A Proofs of Lemmas 15.11 and 15.21 


Lemma f5.il 


Proof. Consider an optimal allocation d = (di,. .., dn), and assume without loss of generality that ui(di) > 
• •• > Vn{dn), and (by the monotonicity of valuations) that = c. Let t = [^1’ N = 

{1,2,..., f}. Without loss of generality, suppose c® — Ylk&N df > ^ n > t. Then, there 

exist ki,... ,kra G Af \ N such that, for all i G [m], d\, > -. We define another allocation 

d = (di,..., dn) G TT^ as follows: let d = (d^ : k G N), and for i G [m], set 

d^ if A: G N, 

0 lik = ki, 

^ • 6®^ ^ if k ^ N and k ^ ki. 

^N,d. ’ 

Note that 

^ dfc < ^ {dl + 6 ^ 3 ) < ^ d{. -I- (n - 3 - d{,, < c® - ^ d\ (37) 

k^N k^N, k^ki k^N k&N 



that is, d is a feasible allocation. Let 

°N,A 

follows using Eqn. (371 that Y^k^N '^1 - 


E'-is 

k^N 


E 

k^N, ky^ki 


dl + hi 


NA 


= (n 


iV,d 


for k G Jf\{N U {ki}) and = 0 for k 

.,2 ^k^N, k^ki (dl + ^N,d^ 2 

~ G - i —- Ji - <{n-t) 

C Yk£N dk 


ki. Then it 


(38) 
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Thus d G 5. Finally note that, for all i G [m], Vk{dki) < \ Yhk&N'^kidk)^ from which we get by the 
monotonicity of the valuations that u(d) = T^keNMdk) + Efc^TV, > J2kMdk) - 

E™i^fc(4J>(l-6)0PT. □ 

Lemma 1131 

Proof. We first observe that, due to the way the valuations are defined in Q, we may assume for fhe purpose 
of compufing an opfimal allocafion d* fhaf V = (J^i, Dk- Indeed, suppose fhaf d* = (d^,..., d*) G j*, 

where d* = (d^ : k £ N), df, 0 D^i for some k' G N, and d\ = (j'E ^JVd* ■ * ^ I™-]) k ^ N. 
Then lef us define a new allocafion d as follows: for each k £ N, we choose dk G Dk such fhaf dk ^ d^ anc 


fhaf d*’ hence u(d) > r;(d*). We nofe furfhermore fhaf d G since for all i, we have 


Vk{dk) = Vk{d\)-, we sef d = (d^ : k £ N), and for k ^ N, define = (r^ • : i £ [m]). Nofe by (14 


Y,{4 - <■•) = 11(4- <ii') + E E - 4 

k'^N ^ kGN 


kGN 


= E<4-4-‘) 1-7^ 


k^N 


Yl,k'^N ^k' 
n- |iV|)2 


< 0 , 


since d^ < d^’* and Ylk'^N'^'k' < for *■ H follows fhaf Yhkdk < Ylkdk* < c, and hence 

d G 5jy 3 as claimed. 

To maximize over S, wifh fhe resfricfion thaf V = |J^ Dk, we iferafe over all subsefs C AA of size af 
mosf — and all parfial selections d = {dk £ Dk : k £ N). For each such choice {N, d), we use dynamic 


programming fo find argmax^g^^ _ r’(d)- Lot d f’o as defined in (14 1 . Wifhouf loss of generalify, assume 
M \ N = {1,..., n — f}. For k £ J\f \ N and r = {r ^,..., r"*) G {0,1,..., (n — |At|)^}”*, define U{k, r) 
fo be fhe maximum value obfained from a subsef of users {1,2,...,A:} C J\f \ N, wifh user j £ [k] having 
demand d* = r®- • for f G [m], where r®- G {0,1,..., (n — |At|)^}, and such fhaf Eje[A:] r"} ^ r-®. 
For fwo vectors x,y £ M®®®, lef us denofe by x * y fhe vector with components (xiyi,..., XmlJm)- Define 
U (1, r) = —oo, if r ^ 0. Then we can use fhe following recurrence fo compufe U{k, r): 


U{l,r) = maxui(6jv j* r) 

r ’ 

f7(A: + l,r) = max {ufc+i(6^ j=t= r^+i) + f7(A:,r - r^+i)} . 

rk+i<r 

Nofe fhaf fhe number of possible choices for r is af mosf n^®®®, and hence fhe fofal fime required by fhe dynamic 
program is Finally, given fhe vector r fhaf maximizes U{n — |N|, r), we can obfain (by fracing back 

the optimal choices in the table) an optimal vector = (r^,..., r^), for each k £ M \ N. From this, we 
get an allocation d G 5, by defining dk = dk for k £ N and, for k ^ N, we choose dk £ Dk such fhaf 
dkfirk* 67V,d and Vkidk) = Vk{rk * 6Ar,d)- □ 
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